Photon Emission from a Medium-Modified Shower Evolution 
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Photons from the interaction of a highly energetic jet with a thermal medium are an important 
contribution to the total photon yield measured in ultrarelativistic heavy-ion collisions and also an 
important probe to study the medium degrees of freedom. Previously this contribution has often 
been computed in the context of a leading parton energy loss approximation. In this work, jet- 
medium interaction photons are instead estimated using a medium-modified shower evolution model, 
where the energy degradation due to vacuum radiation prior to medium formation, the virtuality 
evolution of intermediate states and the photon emission from subleading shower partons is taken 
into account consistently. The results indicate that the leading parton energy loss approximation 
does not appear to work well for photon emission from jet-medium interaction. 

PACS numbers: 25.75.-q,25.75.Gz 
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I. INTRODUCTION 

Photonic measurements constitute an important part 
of the set of observables considered to study the prop- 
erty of the Quantum Chromodynamics (QCD) medium 
created in ultrarelativistic heavy-ion collisions. As with 
dileptons, the other major class of electromagnetic (e.m.) 
observables, their value arises from the relative smallness 
of the e.m. coupling a em as compared to the strong cou- 
pling a s which implies that any photon or lepton has a 
mean free path two orders of magnitude larger than the 
constituents of stronlgy interacting matter. Thus, once 
produced in a medium, e.m. probes escape practically 
without any re-interaction from the medium. 

Thermal photon emission as a tool to study the space- 
time emission of heavy-ion collisions has been studied as 
early as in the SPS era [l|-[4| and has now reached a high 
degree of sophistication, involving state-of-the art mod- 
elling of the medium evolution in terms of event-by-event 
fluctuating or viscous hydrodynamics [5|-|7[ aiming to re- 
produce the measured angular momentum anisotropy co- 
efficient l>2 Q- 

In a different momentum regime, hard direct photons 
produced from perturbative high Pt partonic reactions 
in the medium are often considered as the 'golden chan- 
nel' to measure parton-medium interactions, as measur- 
ing the photon allows to constrain the hard reaction kine- 
matics very well |9l-tl3j. 

In addition to thermal photons originating from the 
medium and hard photons produced in high Pt reactions, 
there arc additional sources of photons related to the in- 
teraction of hard partons with the medium, for instance 
the so-called jet conversion photon contribution [14j or 
the induced radiation contribution [15j . Such photons are 
a relevant background to thermal photon observables, for 
instance they tend to dilute any photon v^ contribution, 
but also dilute the clean kinematics of a 7-triggered hard 



'Electronic address: thorsten.i.renk@jyu.fi 



correlation. However, since jet-thermal photons carry in- 
formation about jet-medium interaction, they are inter- 
esting in their own right and for instance jet -tag ging has 
been suggested to isolate their contribution [16| . 

In [171 ] , the various sources of photons have been com- 
puted and classified as follows: Prompt-direct photons 
directly result from the hard process itself. In contrast, 
prompt fragmentation photons are radiated in the final 
state parton shower of a quark or gluon produced in a 
hard process. Jet-QGP photons come from the interac- 
tion of hard partons with the medium such as the con- 
version reaction, while thermal photons are produced by 
interactions among medium partons. 

This classification is motivated by the energy loss pic- 
ture of the interaction of hard partons with the medium, 
in which the passage of a high Pt parton through the 
medium is treated as energy loss from the leading parton 
due to induced radiation, followed by vacuum fragmenta- 
tion of parton shifted in energy outside the medium (see 
e -g- [18l Il9|). However, a more modern understanding, 
driven by the need to understand fully reconstructed jet 
observables, is to see the energy loss picture an approxi- 
mation to a medium-modified shower evolution in which 
vacuum and medium-induced radiation are not separable 
emission by emission but only on average (see e.g. [201 ] for 
a detailed discussion of the energy loss approximation) . 
There are now several Monte-Carlo (MC) codes available 
to compute in-medium shower evolution |2l|-|26j. 

Since photon emission through Quantum Electrody- 
namics (QED) processes is subject to the same radiation 
phase space considerations as the vacuum or medium- 
induced partonic QCD radiation, this changed view has 
pronounced implications also for jet-medium photons. 
In particular, in a medium-modified shower picture one 
can not separate prompt fragmentation from medium- 
induced radiation photons, and the uncertainty relation 
implies that a significant part of the QCD vacuum ra- 
diation occurs before medium-induced photons can be 
generated, not after the hard parton exits the medium 
as assumed in the energy loss approximation. The aim 
of this work is to explore jet-medium photon emission in 



a medium-modified shower picture, here exemplified by 
the MC code YaJEM [H la. 



II. QUALITATIVE CONSIDERATIONS 

As most in-medium shower codes, YaJEM is based on 
the vacuum shower code PYSHOW 27] which simulates 
the QCD evolution of a highly virtual initial parton as 
a series of 1 — > 2 splittings of a parent parton into two 
daughters with decreased virtuality. The evolution is car- 
ried out until a lower non-perturbative virtuality scale 
0(1 GeV) is reached, at which point a non-perturbative 
hadronization model (such as the Lund model |2g|) is 
used. The vacuum fragmentation photon contribution 
can in this picture be obtained by allowing in addition to 
the QCD processes q — >• qg g —$■ gg and g — ► qq also the 
QED 1 — > 2 splitting q — > qj, i.e. any parton which is 
charged can radiate electromagnetically if there is phase 
space available. 

In a schematic way, jet-medium interaction photons 
can be classified as being produced in two different pro- 
cesses: Bremsstrahlung photons are produced when in- 
teraction with the medium makes a radiation kinemati- 
cally possible by creating radiation phase space (on-shell 
charges can not radiate real photons, but interaction with 
the medium can move the charge off-shell). Diagrammat- 
ically, the simplest such processes are 2 — > 3 reactions in 
which a photon line is attached to a QCD 2 — > 2 scat- 
tering diagram. In contrast, conversion photons are the 
results of the 2 — > 2 processes qg — >■ qj or qq — > gj (with 
qq — > 77 much suppressed due to the weakness of the 
e.m. coupling) in which partons change their identity 
rather than radiation phase space. 

In previous works in the context of the energy loss 
approximation (e.g. [lj, [la, ll3|)i it nas been assumed 
that bremsstrahlung and conversion photons can be com- 
puted, as the medium-induced radiation, for an on-shell 
quark with initially the full energy coming from the hard 
process, and that the fragmentation photon component 
can then be obtained by vacuum-fragmenting the quark 
after its energy has been shifted down due to its interac- 
tions with the medium. 

However, quarks are created with a high initial virtu- 
ality. The Heisenberg uncertainty principle suggests that 
parametrically a virtual state lives for a time r = E/Q 2 
where E is the energy and Q the virtuality of the quark. 
Taking, as assumed by vacuum parton showers, a lower 
perturbative scale of 1 GeV, this implies that even at 
RHIC 20 GeV quark jets have a parametric lifetime of 
~ 4 fm during which the virtuality is above 1 GeV, i.e. 
leading shower partons may be quite far from being on- 
shell while they pass through the medium for a significant 
time. Since initial virtualities parametrically are Q ~ E, 
the first branchings take place at timescales C 0.1 fm, 
i.e. before any medium can form. Thus, typically vacuum 
radiation degrades the energy of the leading quark even 
before it can interact with the medium, and the state 



encountering the medium, rather than being a single on- 
shell quark, is a shower of still highly virtual quarks and 
gluons. In particular, in such a shower there may be sec- 
ondary charges, i.e. also a gluon jet gives in principle rise 
to photon production. 

This has a weak impact on computations of leading 
hadron production, as what matters for the observation 
of the leading hadron is mainly the energy of the leading 
quark before hadronization, and the time ordering (i.e. 
whether vacuum radiation reduced the energy first and 
medium-induced radiation later or vice versa) is irrele- 
vant. This is one of the reasons the energy loss approxi- 
mation works so well for leading hadron production [20j . 
That the same approximation holds is less obvious for 
photon production — it can be argued that for the bal- 
ance of fragmentation photons and bremsstrahlung pho- 
tons a similar argument holds, and that increased energy 
loss implies stronger QCD and stronger QED induced 
radiation, compensated by a fragmentation with reduced 
energy later. However it is not obvious that the momen- 
tum dependence of both contribution cancels. 

A strong impact is however expected for the conversion 
photon channel: Here, the time ordering, i.e. whether 
some energy has been lost before the medium is encoun- 
tered or after is crucial — the rate of conversion photons 
for given momentum is much reduced if the quark en- 
ergy is degraded before the medium forms. Likewise, the 
conversion processes are enhanced by t and u channel sin- 
gularities for near on-shell quarks [l4j which are assumed 
to be screened by thermal masses O(gT) (with g the cou- 
pling constant and T the medium temperature). If the 
singularities are screened by the much larger virtuality of 
intermediate shower states instead, the conversion cross 
section is much reduced. 

In the following, we will investigate the strength of the 
individual contribution in the context of the in-medium 
shower code YaJEM. 



III. THE MODEL 

As described above, the MC code YaJEM is based on 
the PYSHOW algorithm [27| which in turn is part of 
PYTHIA 29]. It simulates the evolution from an initial 
parton with virtuality Qi to a shower of partons at lower 
virtuality in the presence of a medium. In the absence of 
a medium, YaJEM by construction reproduces the results 
of PYSHOW. A detailed description of the model can be 
found in [H, [H, [3(1 . Here the version YaJEM-DE is used 
31] which is one of the best-tested theoretical models 
available for in-medium shower evolution and gives a fair 
account of a large number of high Pt observables both 
at RHIC and LHC [13, EJIH] . 

The medium-modification of the shower evolution is 
implemented via the modification of radiation phase 
space. The medium itself appears only through transport 
coefficients q, e which influence the kinematics (in terms 
of energy E a and virtuality Q a of intermediate virtual 



states a as 
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can be more than one conversion in a shower, or there 
can be both conversion and fragmentation photons from 
the same shower, although in practice this is exceedingly 
rare). 
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IV. RESULTS 
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where r^ is the time at which a fluctuation a is cre- 
ated T a is the lifetime of the virtual state as given by 
a randomized evaluation of the Heisenberg uncertainty 
relation. 

Individual scatterings with the medium are not re- 
solved in this framework. This makes the computation of 
the combined vacuum and bremsstrahlung photon con- 
tribution straightforward by permitting in addition to 
QCD also QED branchings during the shower evolution. 
However, the conversion photon contribution can not be 
computed without resolving the medium. 

In order to get an estimate for the importance of the 
conversion photon contribution, we assume in the fol- 
lowing that for the purpose of evaluating the conversion 
processes only, the medium can be locally decomposed as 
a free gas of quarks and gluons at temperature T (note 
that such an assumption does not lead to a good descrip- 
tion of the measured nuclear suppression factor Raa [35| , 
it is therefore at best an estimate). 

In the following, we utilize the main result for the spec- 
trum of conversion photons, given the distributions f q 
and fq of hard quarks propagating through a medium 
with temperature T given in [14| as 
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Note that the momentum distribution of conversion 
photons is given by the distribution of converting quarks. 
Inserting the distribution f q ( q )(p) — (27r) 3 <5(a; — Xo)5(y — 
yo)S(z — ct)S 3 (p), i.e. evaluating the expression for a sin- 
gle propagating virtual quark state at known position 
and momentum while setting the mass scale m to the 
virtuality m = Q allows to compute the the conversion 
probability of this state per unit time. This can be inte- 
grated from t° to r° + r a to find the probability of the 
conversion during the lifetime of the virtual intermediate 
state, which can be sampled in YaJEM for each virtual 
quark state during the shower evolution. In the case of 
a conversion reaction, the propagating quark is changed 
into a photon and the shower evolution is carried out 
further for the remaining partons (i.e. in principle there 



A. Medium-modified fragmentation photons 

As found in 23], the medium-modified fragmenta- 
tion function computed in YaJEM obeys a scaling law 
that makes the shower largely independent of the pre- 
cise functional form dependence of the transport coeffi- 
cient q(Q on the actual path £ but only on the integral 
AQ 2 = J d(q(Q for most paths which can occur in a 
realistic medium evolution. This allows to characterize 
the strength of the medium modification by the single 
parameter AQ 2 . 
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FIG. 1: Medium-modified fragmentation photon yield from a 
20 GeV parton evolved with YaJEM in vacuum and for two 
in-medium path with different strength of the total parton- 
medium interaction. Solid lines denote quark results, dashed 
lines denote gluon results. 

Fig. Q] shows computed medium-modified fragmen- 
tation photon distribution originating from 20 GeV 
shower- initiating partons for two different strength of the 
medium evolution. Here, AQ 2 = 10 GeV 2 approximates 
a long path through the center of a medium created at 
RHIC conditions, whereas AQ 2 = 5 GeV 2 represents 
the average in-medium path. Overall, the effect of the 
medium is rather modest and largely confined to low E~ l 
which makes it difficult to observe (when convoluted with 
a pQCD parton spectrum as in the computation of a pho- 
ton yield, dominantly the region above 15 GeV where 
differences vanish is probed). 

Fragmenting gluons give rise to a small photon compo- 
nent which comes from the splitting g —> qq early on in 
the shower evolution where one of the daughter quarks 
later undergoes an e.m. emission. The fact that a QED 
splitting is not possible from the shower initiator much 
depletes the distribution at high E 1 . 



In order to bring out the differences between medium- 
modified and vacuum distribution better, we plot their 
ratio (Iaa) in Fig. [2] 
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FIG. 3: Yield of conversion photons resulting from a 20 GeV 
quark for two different paths through the medium, compared 
with the vacuum fragmentation photon yield. 



FIG. 2: Ratio of medium-modified over vacuum fragmenta- 
tion photon yield from a 20 GeV parton evolved with YaJEM 
for two in-medium path with different strength of the total 
parton-medium interaction. Solid symbols denote quark re- 
sults, open denote gluon results. 

Again, the dominance of soft medium-induced e.m. ra- 
diation is clearly visible. Interestingly, the relative en- 
hancement of the medium- modified yield is about a factor 
two higher for gluons, reflecting the fact that gluons al- 
ways radiate photons by splitting into a qq pair. At high 
Pt the relative enhancement is consistent with unity. 

The medium-induced radiation spectrum depends on 
the energy scale of the shower initiating parton. Repeat- 
ing the computation for fragmenting 200 GeV quarks and 
gluons yields no visible in-medium enhancement of the 
photon production within statistical errors. 



B. Conversion photons 

Using the prescription outlined above, the estimate 
for the conversion photon yield from 20 GeV shower- 
initiating quarks is shown in Fig. [3] and compared with 
the vacuum fragmentation yield. 

Since there is no reason to assume that the conversion 
photon yield obeys a scaling law, the estimate is done for 
two different actual paths through a 2+ Id ideal hydrody- 
namical simulation of RHIC central heavy-ion collisions 
— one path from the medium center, and one from the 
periphery (i.e. a vertex displaced by 4 fm towards the 
surface) . 

The conversion photon yield is found to be broadly 
distributed in energy (note that the results of [1J, [l6| or 
the direct evaluation of Eq. ^ for the shower-initiating 
quark would expect it to be of the form ~ 8{E 1 — 
20 GeV)), indicating the strong role of vacuum radia- 
tion prior to the medium interaction. Furthermore, the 
yield is found to be smaller than the fragmentation yield 
throughout the whole energy range. While in [14| a 



strong conversion yield has to be obtained, this yield in 
the investigation here is much reduced in the computa- 
tion here due to the screening of the singularities by a 
large timelike virtuality > 1 GeV for virtual states as 
compared to screening by a thermal mass of order ~ gT. 

Somewhat surprisingly, the yield of conversion photons 
is relatively independent of the actual in-medium path. 
This is the result of an accidential cancellation between 
the T dependence of the conversion probability favour- 
ing a high-density medium and energy degradation by 
medium-induced radiation prior to conversion, favouring 
a lower density medium. 

From Eq. ([3]) one can expect that the shape of the pho- 
ton spectrum reflects the shape of a quark distribution. 
This interpretation is blurred by the fact that the conver- 
sion photon yield is obtained as the result of a spacetime 
integration over an evolving medium, however the result- 
ing shape is found to compare well with the quark distri- 
bution in the shower after an evolution time of about 2 
fm. This timescale agrees reasonably well with the peak 
strength of medium-induced radiation (earlier times are 
screened by interference, later times are suppressed due 
to the rapidly decreasing medium density). 



V. DISCUSSION 

The results obtained above suggest that the energy loss 
approximation can not be expected to give sufficiently 
precise results for photon emission from jet-medium in- 
teraction. In particular, the broad distribution of con- 
version photons seen in Fig. [3J given a quark with a well- 
defined initial energy suggests that jet-tagging of conver- 
sion photons as suggested in [16( might not be feasible, 
and that the relative strength of the conversion photon 
contribution to the total should be re-evaluated in models 
beyond the energy loss approximation. 

Fig. Q] suggests that overall the strength of the medium 
modification of photon emission is not large, and that the 



vacuum fragmentation contribution alone is a reasonably 
proxy for the total. A detailed investigation of the precise 
phenomenological consequences of this work is left for a 
future investigation. 
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